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Abstract 

We discuss form factors and coupling constants for the 7*7r°7, pint and ^irp interactions gener- 
i— l ated by a model field theory that produces finite size qq meson modes. The approach implements 

dressing of the vertices and propagators consistent with dynamical chiral symmetry breaking, gauge 
invariance, quark confinement and perturbative QCD. 



^ 1 Introduction 

Q\ We generate mesons as qq composites with intrinsic size from the Global Color Model (GCM). [1, 2] 
^ This is a phenomenologically successful QCD-based model field theory that permits an accessible and 



If} covariant description of meson substructure and interactions in terms of dynamically dressed quark 

o 

0\ 



^ degrees of freedom. The GCM action is 



S[q, q} = J d A x q(x)( 7 ■ d x + m)q(x) + ± | d'xd'y j^D^x - y)f v {y) (1) 

^ where j^(x) = q(x)^--y tl q(x). The formulation is in Euclidean space with metric 6^, and in this work, 
g we consider only two flavors [u and d). The GCM represents the gluon sector by a finite range effective 
gluon two-point function (D^(x) = S^D^x), in the Feynman-like gauge chosen here), thus formalizing 
the Abelian approximation to QCD. Nevertheless, the chiral anomalies [2] are properly embedded and 
this is important for several of the topics addressed here. A closely related model is that of Nambu and 
Jona-Lasinio [3] where the interaction is a contact one. In practical terms the difference is important. 
Quark confinement can be incorporated within the GCM and will prevent a spurious qq width for the 
heavier mesons. Finite size effects associated with dynamical quark self-energies and meson Bethe- 
Salpeter (BS) amplitudes provide natural regularization of many of the loop integrals for key physical 
quantities. 

Studies of mesons and baryons within the GCM require hadronization techniques [4] and a recent 
review [5] is available. For mesons, a brief summary is the following. A functional change of path 
integration variables allows the second term of (1) to be eliminated in favor of q(x)A a B a (x, y)q(y) and 
a term quadratic in B. These auxiliary boson fields B a (x,y) transform as q(y)A a q(x) where A a are 
the matrices from Fierz reordering of the current-current term. The saddle point configuration of B a 
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yields a quark self-energy (equivalent to the ladder Dyson-Schwinger equation (DSE)) and expansion 
about there identifies a meson action in terms of field variables B for propagating modes. The resulting 
bare inverse propagator may be used to generate eigenvectors T(q, P) for expansion of the bilocal fields 
as B(q,P) — > T{,(q, P)b(P), leaving an effective local field variable b(P) for each meson mode. On the 
mass-shell, the T so defined are solutions of the ladder Bethe-Salpeter equation. The momentum q is 
conjugate to x — y and P is conjugate to 

To allow for electromagnetic (EM) coupling to the meson modes, the bosonization method can be 
generalized [6] to account for a background electromagnetic field minimally coupled (d^ — > d v — iQA„) 
to the bare quarks. The saddle point configuration jf?o[^] now dresses the photon-quark vertex as well 
as the quarks. Gauge invariance of the EM coupling at the quark level can be translated to the meson 
level. [6] The method may be adapted for ease of computation, and for our purposes here, the following 
simplified version of the resulting meson action shall suffice 

oo / \ n 

S[A,ir,p,...] = Tr£^^[£(i 75 f -Trr^ + ^f -^r p + "Of 

+ l I ll^D-^ir + ptpD^tpp + •••] + 6S[A,p,u]. (2) 

The trace here is over spin, flavor and color as well as space-time coordinates. We have employed only the 
dominant covariant (i.e. single Dirac matrix) for each meson mode. The EM field A v appears only in the 
dressed quark propagator •S'fAJ and in the term SS[A, p, cj] which is linear in the indicated meson fields 
and is at most linear in A v in the more useful approaches. When A v = 0, the n = 1 contribution from 
the first term cancels the SS term since in that case the meson fields are defined as fluctuations about 
the saddle point configurations. Then S = S[0] has ladder DSE content giving S(p)~ 1 = i~f -p+S(p)+m. 
When A v ^ 0, rather than expand about the new saddle point configurations, we choose to expand about 
the configurations that produce S(p, &) _1 = S(p — k)S(p)~ 1 +r„( £ |^;p — k)A u {jp — k) where T„(jp\q) is 
a dressed quark-photon vertex that is convenient for calculations. The third term of Eq. (2) depends 
on the choice of Y u but it does not contribute to the interactions we study here. 

A central element of the resulting effective action for bare hadrons is the dressed quark propagator. 
The necessary calculation of meson propagators and interaction vertex functions is complicated by the 
fact that the dressed quark propagator must be continued to complex Euclidean momenta to access the 
meson mass-shells in the timelike domain. In order to explore the capabilities and limitations of the 
overall approach, recent efforts have chosen the dressed quark propagator as the vehicle that carries 
the phenomenological input. Within the representation S(p) = —i-y -p crv(p 2 ) + a s(p 2 )i we employ the 
parametrized amplitudes [7] 

*s{x) = c + , 6 , 6 {bo + + -^(1 - e- 2 ^)) (3) 

bix b 3 x v b 4 x 1 x + m l 

. 2(x + m 2 ) - 1 + e - 2 (-+™ 2 ) _ 2x 

= ~ — 2(xT^y cme ' (4) 

with x = p 2 /A 2 , as = Acr^, a v = A 2 oy, m = m/A, where m is the bare quark mass and A is the 
momentum scale. This S(p) is an entire function in the complex momentum plane, a condition sufficient 
to ensure the absence of quark production thresholds in S-matrix amplitudes for physical processes. [8] 
This form is guided by a confining model DSE [9] and by the behavior found in realistic DSE studies. [10] 
It is consistent with pQCD in the deep Euclidean region up to ln(p 2 ) corrections. The parameters are 
A = 0.568 GeV, c = 0.0406, m = 6.7 MeV, and (b ,b u b 2 ,b 3 ,b 4 ) = (0.118,2.51,0.525,0.169,1 x 10" 4 ). 
The fitted soft chiral physics quantities produced by this parameterization are [7] /„- = 92.4 MeV, 
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< m >iGeV 2 = (232 MeV) 3 , = 140 MeV, r e ™ = 0.54 fm together with reproduction of the 

experimental tttt scattering lengths to within 20%. 

The choice of quark propagator parameters is equivalent to an implicit choice of effective gluon 
propagator and vertex in the underlying DSE dynamics. [10] Ghost contributions in DSE studies have 
been studied in Landau gauge and shown not to modify the qualitative features of quark and gluon 
propagators. [11] Quantitatively, ghosts provide a small (< 10%) effect. We therefore expect that a well- 
chosen phenomenological quark propagator captures those aspects of the underlying dynamics that is 
of practical importance for the chiral physics quantities that guide it. Independent applications, such 
as those contained here, provide important additional tests of this standpoint. 

In terms of the equivalent representation S(p)~ 1 = i~f • pA(p 2 ) + B(p 2 ) + m, the chiral limit for the 
pion mass-shell Bethe-Salpeter amplitude is Y v {jp]P 2 = 0) = i^^rB{jp 2 ,m = 0)//„- since the DSE for 
B(p 2 ) and the Bethe-Salpeter equation for the pion invariant amplitude become identical. [12] Hence 
the pion is both a qq bound state and a Goldstone boson. For finite m we use the pseudoscalar form 
r„-(p;P) i^^rB{jp 2 ,m)/ /„- consistent with PCAC. We use the Ball-Chiu [13] Ansatz for the dressed 
quark photon vertex which is T^pj^) = Q^ii(p;q), where Q = |(r 3 + |) is the quark charge operator, 
and 

Up;?) = -i^\{A( P+ ) + A(p_)) + -^-[i 7 • p(A(p_) - A( P+ )) + (B(p_) - B(p + ))} (5) 

with p± = p ± |. This vertex satisfies the Ward-Takahashi identity, transforms correctly and has the 
correct perturbative limit; but it is not unique. We have investigated several other choices [10] and 
found no significant change in our conclusions. 

2 The 7r°77 Form Factor 

The pion charge form factor for space-like momenta is one of the simplest but non-trivial testing 
grounds for applications of QCD to hadronic properties. A closely related quantity that has received 
less attention is the 7*7r° — > 7 transition form factor. [14] Here the photon momentum dependence maps 
out a particular off-shell extension of the axial anomaly. [15] Presently available data for this transition 
form factor in the space-like region Q 2 < 2.5 GeV 2 is from the CELLO [16] collaboration at the PETRA 
storage ring where the process e + e~ — > e + e~7r° was measured with geometry requiring one of the two 
intermediate photons to be almost real. There is renewed interest in this transition form factor due to 
the prospect of obtaining higher precision data over a broader momentum range via virtual Compton 
scattering from a proton target at CEBAF. [17] Features distinguishing the present work [18] from a 
previous quark loop study [19] are: quark confinement (thus eliminating spurious quark production 
thresholds), dressing of the photon-quark vertex, and the dynamical relation between the pion Bethe- 
Salpeter amplitude and the quark propagator. The latter elements are crucial for producing the correct 
mass-shell axial anomaly independently of model details. 

The relevant part of the action (2) is identified by expansion to second order in A v . Only the first 
term (with n = 1) contributes and we find t>[7r°77] = — Tli[SY ^A^SY U A U SY ^-k ]. That is 

r d 4 Pd 4 

S[tt 77] = / -(^-M-P - Q)A,(Qy(P)A^P, Q), (6) 
where the vertex function is given by the integral 

A„(P, Q) = -tr j ^ S(k-P- Q)V.{k - I - |; -P - Q)S(k) 

xY u {k - |; Q)S(k - Q)T«(k ~\-Q\ P). (7) 



3 



Figure 1: The quark triangle diagram for the generalized impulse approximation to the 7*7r°7 vertex. 

The momentum assignments are shown in the quark triangle diagram of Fig. 1. The general form of 
the vertex allowed by CPT is A^(P,Q) = —i^j-e^apPaQ p g(Q 2 ,P 2 ,P • Q) where e 4 i 2 3 = 1, a. is the 
fine-structure constant, /„- is the pion decay constant, and g is the off-mass-shell invariant amplitude. 
With the one photon mass-shell condition (P + Q) 2 = 0, the invariant amplitude, denoted by g(Q 2 , P 2 ), 
is the object of the present work. For a physical pion the shape of the 7*7r° — > 7 transition form factor 
is given by g(Q 2 , — m 2 ). The chiral limit for the physical 7r° — > 77 decay amplitude is fixed at -^j- by 
the axial anomaly [20] which gives an excellent account of the 7.7 eV width and requires (7(0,0) = 1/2. 
This follows only from gauge invariance and chiral symmetry in quantum field theory and provides a 
stringent check upon model calculations. 

The 7r°77 vertex function in (7) is now completely specified in terms of the quark propagator. An 
exactly parallel situation holds for the spacelike pion charge form factor and this model has previously 
been shown to provide an excellent description of the data. [7] No adjustment of parameters is made 
in the present application. At Q 2 = 0, our numerical evaluation yields g^o^ = g(0, —ttl 2 ^) = 0.496, in 
agreement with the previous application of this model [7] and in good agreement with the experimental 
value 0.504 ± 0.019. The chiral limit of this approach has been shown [7] to correctly incorporate the 
exact result (7(0,0) = 1/2 produced by the axial anomaly independent of the form and details of the 
quark propagator. 

The obtained form factor F(Q 2 ) = g(Q 2 , — m^.)/^(0, — m 2 ) at the pion mass-shell is displayed as 
Q 2 F(Q 2 ) in Fig. 2 by the solid line along with the CELLO data. [16] We calculate a "radius" or 
interaction size, defined via (r^ 7r o 7 ) = —6 F'(Q 2 )\q 2=0 , of 0.47 fm while a monopole fit to the data 
yields [16] 0.65 ± 0.03 fm. The insert compares our low Q 2 result with a recent result from a QCD sum 
rule approach [21], and a monopole form [22] that interpolates from the leading asymptotic behavior 
F(Q 2 ) — > 8TT 2 f 2 /Q 2 argued from the pQCD factorization approach. [23] In the latter two approaches 
there is ambiguity due to: A) the unknown momentum scale at which perturbative behavior should 
set in; and B) assumptions for the pion wavefunction and how it should evolve with the momentum 
scale. [21] Within the present approach, both the photon coupling and the pion wavefunction evolve 
with Q 2 in a way determined by the evolution of the dressed quark propagators. This produces, in 
a single expression, both the ultra-violet behavior required by pQCD and the infra-red limit dictated 
by the axial anomaly. We note that the employed pion Bethe-Salpeter amplitude B{jp 2 , m)/ /„- has 
the correct leading power law behavior mA 2 /p 2 /„- which implements the hard gluon contribution that 
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Figure 2: The 7*7r°7 transition form factor. 

dominates pQCD. 

The dotted straight line in Fig. 2 is the pQCD factorization [23] limit Q 2 F(Q 2 ) -> 8tt 2 / 2 = 0.67 GeV 2 . 
Although our non-factorized calculation reaches this value near Q 2 = 3 GeV 2 , we find a slow decrease 
with higher Q 2 consistent with a logarithmic correction. An excellent fit to the numerical results for 
3.3 GeV 2 < Q 2 < 10 GeV 2 is provided by F(Q 2 ) = A [1.0 + B Q 2 ln(C Q 2 )]" 1 , where A = 1.021, B = 
0.461 /m 2 = 0.777 GeV" 2 and C = 1.16 /m 2 = 1.45 GeV" 2 . (We have neglected the anomalous dimen- 
sion of the quark propagator in this calculation, which would modify the power of the ln-correction.) 

The logarithmic correction to the anticipated 1/Q 2 asymptotic behavior can be attributed to the 
persistent nonperturbative nature of the coupling to the final state soft photon in this exclusive pro- 
cess. [18] Numerically we find that, if a bare coupling were to be used for both photons as is implicit 
in the pQCD factorization approach, F(Q 2 ) would eventually approach 8TT 2 f 2 /Q 2 . The turn-over in 
Q 2 F(Q 2 ) near 3 GeV 2 predicted in Fig. 2 is barely within the Q 2 limit of 4 GeV 2 anticipated for 
measurements at CEBAF if a 6 GeV electron beam becomes available. [17] This turn-over and the 
logarithmic corrections generated by the loop integral are features also found in the parallel approach 
to the pion charge form factor. [7] For a study of the behavior of the vertex function off the pion mass 
shell, as needed for the anticipated CEBAF experiment, see Ref. [18]. 

3 The pirn Form Factor 

From (2) we identify S [piviv] = — Tr ^^^r • p^T p (Si , -/5T ■ 7?r 7r ) 2 j , which yields 
where the vertex is 

MP ' Q) = I ^Y fp{h + y ; Q)T * {k + ? ; _p " ~ ? ; p " l )T " ( *' p ' Q) (9) 
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with 



T^k,P,Q) = tr 



(10) 



From symmetry properties, it is not difficult to show that A P (P, Q) = — A /J (— P, = A P (P, — which 
requires the general form 

A,(P, Q) = -P,P p7r7r (P 2 , Q\ (P • Q) 2 ) - Q.P ■ QH P ^(P 2 , Q 2 ,(P- Qf). (11) 

With both pions on the mass-shell, (P - f ) 2 = (P + f ) 2 = -m 2 . Equivalently, P ■ Q = and 
P 2 = — m 2 — ^j- so that only the first term of (11) survives. The corresponding form factor Fp^^Q 2 ) 
contains the coupling constant as its mass-shell value, i.e. Fp^^Q 2 = —m 2 ) = g p ^. If the form factor 
is held at this value for all momenta, one obtains the point coupling limit which can be expressed in 
the standard form 

S [p-ww] = -0 pmr J d^xp^x) ■ tt(x) X d ll Tr(x). (12) 





Q 2 (GeV 2 ) 



Figure 3: Behavior of the p-ww form factor with p momentum in the viscinity of the mass-shell at 
-0.6 GeV 2 . 

The previous study [24] within the GCM made the approximation g p7rK ~ P p7r7r (P = Q = 0) in order 
to avoid the occurrence of complex momenta in the arguments of the propagators and vertex functions in 
the integral (9). Complex momenta are unavoidable when a Euclidean formulation is continued to handle 
timelike external momenta. However when propagators and BS amplitudes are produced numerically 
it is a difficult and lengthy matter to extend such solutions to the complex plane. The utility of the 
parameterized representation, (3) and (4), of the quark propagator is that the continuations required 
to properly implement mass-shell constraints as needed for extraction of g p7rK and other quantities are 
a simple matter. Since the propagator parameters have so far been constrained only by pion physics, 
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there is no guarantee that continuation of S(p) away from the real p 2 axis by an amount proportional 
to m 2 will be correct. The present exploratory calculations involving vector mesons constitute a first 
investigation of this matter. 

The investigation here is limited since we have not made an independent calculation of the p BS 
amplitude T p (p). This requires knowledge of the effective gluon propagator D(r) compatible, in a ladder 
DSE sense, with the given quark propagator. Although some progress has recently been made in this 
direction, [25, 26] in the present work we directly parameterize T p (p 2 ) oc e~ p l a . This has proved to be 
effective in a related study. [27] The strength is set by the canonical normalization condition [20] which 
is equivalent to ensuring that the residue at the pole of the p propagator is unity. Only the quadratic 
p dependence from the first term of (2), i.e. the quark loop contribution to the p inverse propagator, is 
needed to implement this. The range a is then adjusted to reproduce the empirical value g e p ^ = 6.2. 
This requires a = 0.208 GeV 2 . The corresponding p — > tvtv decay width, given by 

2 

"^"1712 

is 151 MeV. The calculated form factor Fp^^Q 2 ) is shown in Fig. 3 for timelike and spacelike momenta 
in the viscinity of the mass-shell. The main conclusion from this calculation is that the previous 
approximation of using zero momentum to extract a coupling constant [24] can underestimate the value 
by almost a factor of 2. 



4m 2 



(13) 



4 The 77rp Form Factor 

The empirical p BS amplitude set by g p7rK as above enables a parameter-free prediction for the -ynp 
vertex. Apart from being a consistency check in this manner, the -ynp interaction together with the 
77T7T, 777T, and the p-ww processes, provide important guidance for extending the present approach 
to nonperturbative QCD modeling of meson physics to phenomena not dictated by chiral symmetry. 
Within nuclear physics, the associated isoscalar 7*71-/0 meson-exchange current contributes significantly 
to electron scattering from light nuclei. In particular, our understanding of the deuteron EM structure 
functions for Q 2 2 — 4 GeV 2 is presently hindered by uncertainties in the behavior of this form 
factor. [28]. 

Expansion of (2) to first order in the EM field yields the 771-/0 interaction as the pair of contributions 
S [771-/0] = -TT[ST IJ A IJ Si'y 5 T ■ tvT^Si^t ■ p u T p ] - Tr[5T ^A^Si^f ■ p u Y p Si^ 5 f ■ 7rTV]. (14) 
With a vertex function A^(P, Q) defined by 

S[l7Vp] = ~I d J^f MQ)ji{ ~ P ~ f } " MP ~ f } V(P ' g) ' (15) 
one may combine the two terms in (14) to obtain the integral 

X tr[5(fc + -|)r i ,(fc + ;<J)S(* + + |)i 75 5(fc_)i 7 „]. (16) 

Here k± = k ± y . The jivp EM current ( — (2iv) 4 g J* s ^ ) will be conserved as a consequence of maintaining 
EM gauge invariance at the quark level. Explicitly, use of the Ward-Takahashi identity Q v T v {k; Q) = 
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QS- x (k - f ) - QS- x (k + f ), in (16) immediately gives Q^K^P^Q) = 0. The general form of the 
vertex function can be shown by symmetries to be A^(P,Q) = —i^-e^apPaQ p 5 p7 r 7 f(Q 2 ,P 2 ,P • Q), 
as is expected for a coupling arising from the chiral anomaly. We have used the standard definition of 
the coupling constant so that, at the triple mass-shell point, the form factor / = 1. 

Our numerical evaluation yields g lvp = 0.5. The experimental p + — > 7r + 7 partial width (67 ± 7 keV) 
determines the empirical value = 0.54 ± 0.03. The -ynp form factor obtained with on-mass-shell 7r 
and p, and weighted by Q 2 , is shown by the solid curve in Figure 4. Also shown for comparison is the 
vector dominance model (VDM) phenomenology initially used to include such meson exchange effects 
in electron scattering analysis, [29] as well as the result from a free constituent quark loop with bare 
photon coupling. [28, 30] The quark-based results produce a much softer form factor than is produced by 
the VDM assumption. Above 50 fm~ 2 , which is readily accessible in electron scattering, the differences 
are serious. 
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Figure 4: The photon momentum dependence of the -ynp form factor from the on-shell approximation 
for both mesons. 



5 Discussion 

The photon-quark vertex employed in this work does not contain an explicit vector meson pole. The 
absence of the vector meson dominance mechanism in the explicit photon vertex is, in fact, consistent 
with the tree-level nature of the action (2). The n = 1 part of the first term, in combination with the 
third term, produces direct coupling of the photon with vector mesons. When the analysis is carried 
to first-order in meson loops, EM coupling mediated by a vector meson propagator will be produced 
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consistently as an additional mechanism to that considered here. More detailed exploration of this point 
is the subject of a future article. [31] 

At the photon mass-shell point, the coupling constants obtained here at tree-level for both 7^7 
and jTvp are within 10% of experiment. Our calculation for jivp is equivalent to the th term in the 
meson loop expansion. An estimate of the meson-loop corrections to our result would be interesting. 
At this stage we can only remark that, within the same approach to the pion EM form factor, pion- 
loop corrections were seen to contribute at the level of < 15% to the charge radius. [32] Similarly, the 
pion-loop contribution to the p mass has the correct sign and magnitude to generate most of the p — us 
mass splitting [33] and is a 2% effect. The emerging picture is that a representation of low-mass mesons 
in terms of dynamically dressed quarks may capture the dominant quantum loop effects. Subsequent 
meson loop dressing would have to overcome distributed coupling produced by finite size effects. 

Instead of considering the underpinnings of the present investigations to be a result of bosonization 
of the GCM action to the bare meson level, there is a more general viewpoint. [10] Selected truncation 
of the tower of coupled Dyson-Schwinger equations of QCD, together with use of a generalized impulse 
approximation, points to the same end result for the processes we consider here. However, the road to 
higher order effects is different in each. Progression through the loop expansion in the effective action 
approach defines a rigid ordering of physics content. There is potentially more freedom to develop an 
efficient ordering with the second approach as one moves beyond the impulse approximation and with 
less severe DSE truncation. We have taken the former viewpoint to emphasize the immediate relevance 
to effective hadronic field theory models. It is becoming quite feasible to generate such models with 
many of the previously phenomenological coupling form factors now given a quark basis that, although 
approximate, captures the dominant influence on dynamics at the hadron size scale. 
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